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Abstract 

Based on previously computed parameters for the electron-phonon 
couplings and the Coulomb exchange, we compute and classify the static 
Jahn- Teller distortions, i.e. the minima of the lowest adiabatic potential 
energy surface, of Cg^, for all values of charge 1 < n < 9 and spin. We 
compute the intra-band electronic excitation energies in the different 
optimal geometries in the sudden approximation, and find a spread of 
the electronic states of roughly 1 eV. We also obtain the leading vibronic 
quantum corrections to the ground-state energy, equal to zero-point 
energy lowering due to the softening of the phonons at the adiabatic 
Jahn- Teller minima: these non-adiabatic corrections are so large that 
for 4 < n < 6 states of different spin symmetry turn lower than the 
high-spin adiabatic ground state. 
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1 Introduction 



Low-spin states are associated to larger distortions, thus larger energy gains, 
than high-spin states in degenerate electron-phonon coupled molecules and 
impurity centers. Electron-electron Coulomb repulsion opposes this tendency, 
favoring high-spin states instead, in accord to the first of Hund's rules. The 
Jahn- Teller ( JT) systems Cqq" are no exception to this rule: if electron-phonon 
coupling was the only relevant interaction, then the n-holes ground state would 
be either of spin S — (even n) or S = ~ (odd n). As was recently shown pQ, 
in positive fullerene ions the size of Coulomb interaction is sufficiently large to 
enforce Hund's rule: the ground states of Cg^ was calculated to always be high 
spin [S — | for n < 5, S = for n > 5) in the adiabatic approximation. 
This result is confirmed for n = 2 by NMR investigation of solid-state com- 
pounds j2j. The JT distortions in Cgo~, though strongly counteracted by the 
larger electron-electron repulsion, yet represent an important, and still largely 
unexplored, contribution to the energetics of Cg^ . Investigation of this contri- 
bution, and in particular of the corrections to the adiabatic approximation, is 
the main subject of this work. 

The JT model relevant for Cgg" is conventionally indicated as h n <g> (A + 
G + H), where h refers to the fivefold-degenerate highest occupied molecular 
orbital (HOMO), and A, G, H refer to the 2 nondegenerate A g , 6 fourfold- 
degenerate G g and 8 fivefold-degenerate H g molecular vibration modes that are 
linearly coupled to the h u states according to icosahedral symmetry [TJEJIl]. 
We investigate this model by treating the normal coordinates for these vibra- 
tional modes as classical variables, and searching the minima of the adiabatic 
potential energy surface in the 66-fold dimensional space of these distortions. 
Each of these static JT configurations is characterized by a reduced symmetry 
from icosahedral to some (usually) lower symmetry. New vibrational frequen- 
cies arise at these local minima: we determine these frequencies by evaluation 
of the Hessian matrix at the minimum The lowering of the vibrational 
frequencies gives the leading quantum correction to the adiabatic approxima- 
tion. The original icosahedral symmetry of the problem is restored once the 
presence of several equivalent optimal distortions is recognized, and quantum 
tunneling between these wells is allowed. Proper accounting of tunneling gives 
the next-order quantum correction, but in the present work, we limit ourselves 
to the study of the local properties of the wells and the connectivity of the sets 
of minima in distortion space, for all values of charge n and spin S. 

The competing intra-molecular exchange of Coulomb origin and the JT in- 
teraction both contribute to the computed spectrum of excitations. Differences 
in energies of the fully relaxed configurations at different spin compare directly 



with spin gaps as could be measured in "slow" spectroscopies such as electron 
or nuclear magnetic resonance. In contrast, the electronic excitation energies 
computed keeping the molecular geometry fixed in the lowest minimum com- 
pare directly with the vertical excitations probed by fast optical spectroscopies. 
Both these class of quantities are reported in this work. 

This paper is organized as follows: Sec. El introduces the model and the 
parameters used in this calculation, which is then described in Sec. El along 
with the properties of the JT minima for all values n and S; Sect. 0] contains 
the vertical excitation spectra. The the zero-point non-adiabatic corrections 
are described in Sec. El The results are discussed in Sec. El and connectivity 
matrices are collected in an Appendix. 

2 The model Hamiltonian 

We report here for completeness the model Hamiltonian previously introduced 
in Ref. pQ to describe the physics of the holes in the h u HOMO of Ceo fullerene: 



where 

H e -v 

H e -e 

are respectively the single-particle Hamiltonian, the vibron contribution (rep- 
resenting the phonon kinetic energy plus the restoring potential expanded 
to quadratic order around the equilibrium configuration of neutral C 60 ), the 
electron- vibron coupling (in the linear JT approximation) [3J El; and finally 
the mutual Coulomb repulsion between the electrons. The c\ m denote the cre- 
ation operators of a hole in the HOMO, described by the single-particle wave 
function tp ma (r). a indicates the spin projection; to and n label the compo- 
nent within the fivefold degenerate electronic HOMO multiplet, based on the 
C 5 quantum number to from the Ih D D 5 D C 5 group chain jU [7] . % counts 
the phonon modes of symmetry A (2 A g , 6 G g and 8 H g modes). C^f are 
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Table 1: Computed mode eigenfrequencies and e-v linear coupling parameters 
of the h u HOMO in Ceo The classical single-mode JT stabilization energies 
E s are tabulated for both D5 and D% distortions, for one hole in the HOMO. 

Clebsch-Gordan coefficients [7] of the icosahedral group Ih, for coupling two 
h u states to phonons of symmetry A. r is a multiplicity label, relevant for 
modes of A = H g symmetry only [HI Ej • QiA/i are the molecular normal- mode 
vibration coordinates (measured from the adiabatic equilibrium configuration 
of C 60 ), and the corresponding conjugate momenta. Spin-orbit is exceed- 
ingly small in C 60 [H] and it is therefore neglected. 

The electron-vibron (e-v) couplings g'- A are conveniently expressed in units 
of the corresponding harmonic vibron quantum of energy hu^. In this cal- 
culation we adopt the numerical values of the e-v coupling parameters, listed 
in Table ^ from the Density Functional (DF) calculation of Ref . [0] , and a 
second calculation P yields couplings in substantial accord with those of Ta- 
ble The numerical factors k Ag = 55, k° 3 = (f) 2 , k Hg = 1 in H c - V have 
been introduced for compatibility with the normalization of Ref. jH]- 
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Table 2: The Coulomb parameters for Cgg", as obtained from the DF cal- 
culations of Ref. pp. One of the tabulated parameters (e.g. F\) is a linear 
combination of the five others. 



The Coulomb matrix elements are defined by: 
w a ^{m,m';n,n') = d 3 r /dV <f* ma (r) ^4v(r') u a y(r, r) <p na (r) yw(r') (6) 



where u a y(r, r') is an effective Coulomb repulsion, screened by the other elec- 
trons of the molecule. Detailed symmetry analysis shows pQ that, assuming 
spin-independence of the orbitals, this set of coefficients can be expressed as 



w a;a ,(m,m>;n,n>) = £ F^ fe<^c£#) 

r,r'A V A 4 / 



(7) 



in terms of a minimal set of independent parameters F r ' r '' A . A DF calculation 
of these parameters was carried out in Ref. pQ , and for our calculation we adopt 
those values of the Coulomb parameters, which we report for completeness in 
Table 121 For the Coulomb parameters we use the shorthands 

F x = F A °, F 2 = F G \ F 3 = F 1,1,Hb , F 4 = F 2 ' 2 '^, F 5 = F 1,2,Hg , (8) 

and the combination 

V 5 45 9 9 J K J 

U defines an average Coulomb repulsion within the n-holes multiplets, so that 

E^{n) = Tr n (H + H c . c ) = en + ~ 1} . (10) 

It should be noted that U differs from the usual definition of the Hubbard U, 
involving the lowest multiplet in each n-configuration: JJ mm = E mm (n + 1) + 
E mm (n — 1) — 2E mm (n). This second definition is inconvenient here, since it 
depends wildly on n. 
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Table 3: The total adiabatic energy ^ ad (Q m i n ) (in meV) of the lowest elec- 
tronic state for each n and S, including the e-e and e-v contributions from 
H vlh + H c -v + # e -e (but excluding the [Un(n — l)/2] term), for Cqq. The last 
two columns distinguish the vibrational (H y ib) and electronic {H e -v + H e - e ) 
contributions. 

3 The adiabatic calculations 

We approximate the vibron operators 0iA/i with classical coordinates, in the 
spirit of the adiabatic approximation. In an orbitally degenerate situation (as 
for the Cqq ions at hand) the adiabatic approximation usually yields fairly 
accurate energetics in the limit of large the e-v couplings, so that tunneling 
between equivalent minima can safely be neglected ^Oj- The phonon kinetic 
term in Q is neglected in the adiabatic approximation. In Sect. |S] we will 
partly restore this term by taking into account quantum zero-point energies. 
In any classical statically JT-distorted configuration, the icosahedral symmetry 
is broken: therefore states of different icosahedral symmetry representations 
are inter-mixed. Only the total number of holes n, total spin S and its pro- 
jection S z are conserved upon distortion. Here, we neglect any change of the 
Coulomb Hamiltonian upon distortion, and we assume therefore that H e - e is 
still determined according to Eqs. (j5|6|7j) by the same parameters Fj of Ta- 
ble El as in icosahedral symmetry. Also, we assume no change of the phonon 
frequencies Ui\ and couplings g r iK upon charging. 

For each n, S and Ms, we allow the 64 (6 x 4 G g plus 8x5 H g ) phonon co- 
ordinates to relax, and determine the optimal distortion, by full minimization 



of the lowest adiabatic potential sheet V^ ad (Q) in the space of all the phonons 
coordinates Q. We leave the A g modes out, since they contribute a trivial 

E As (n) = -\n 2 Y.9LM^ = " L79 meV ' ( n ) 

spin- and symmetry-independent term to the energetics. Because of particle- 
hole symmetry, charges n > 5 can always be reduced to the computed charges 
n < 5. In Table El we report the resulting optimally-distorted energy in 
each spin sector, based on the electron-electron (e-e) and e-v couplings of Cg^ 
ions, as previously published in Ref . pQ . The main outcome of the adiabatic 
calculation is that positive Cgo ions favor high-spin ground states (contrary to 
the analogous finding for negative ions). 

In the present contribution, we extend the previous calculation to obtain 
the complete set of all the equivalent minima for each (n, S) sector. To this 
purpose, we generate about a hundred randomly distributed distortions away 
from the 1^ high-symmetry point, and let the molecule relax to the closest min- 
imum, by combined standard (simplex and conjugate-gradients) minimization 
algorithms. We then apply the symmetry operations of the icosahedral group 
to the each of the minima found, in order to locate any possibly missing min- 
imum. Although the method employed is not deterministic, the symmetry 
analysis makes the probability that any set of minima is incomplete utterly 
negligible. Thus, for each n and S we obtain a complete set of equivalent 
global minima. In the few cases where the minimization leads to non-global 
minima, we have discarded them based on simple comparison of the adiabatic 
energies. 

In Table |U we summarize some global properties of the obtained JT minima 
for all charge and spin states. In these multi-mode JT systems, the local 
symmetry of an optimal distortion is described in terms of the subgroup Gi oca i 
of symmetry operations which leave that minimum invariant. We remind that 
the minima in the n = 1 S — ~ case, where e-e interaction is unimportant, 
were found to be 6, of local D 5 d symmetry [31 HUE]. For 2 < n < 8, where the 
role of e-e interaction is crucial, the number of JT minima follows from the 
local symmetry: it is generally given by the ratio |/h|/|Gio Ca i| of the orders of 
the icosahedral group (120) and of the invariant subgroup. 

Special care has to be taken for n = 5 holes. Here, in addition to the icosa- 
hedral symmetry, the system is particle-hole symmetric, i.e. invariant under 
exchange of fermion creation and annihilation operators. This transformation 
leaves the Coulomb Hamiltonian if e -e invariant, while the the vibron interac- 
tion H e - V is unchanged provided that a sign change of the vibron coordinates 
Qi t A,/j, is also performed. Hence, for a given minimum Q m i n , also its opposite 
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Table 4: The number and the local symmetries of the JT minima for given 
charge n and spin S. In the 5th column the number of neighbors of all orders 
are listed for a given minimum. The last column gives the total amount of 
dimensionless JT distortion at each minimum. 



—Qmin is a minimum of the potential energy surface. In the case n = 5 S = |, 
this leads to a doubling of the minima: the local Civ symmetry would lead to 
30 minima, but 30 more equivalent minima are added in the opposite positions 
by particle-hole symmetry. For n — 5 5" = | instead, the number of minima 
remains 30, since for each minimum there is one of the Ih symmetry operations 
that transforms this minimum into its opposite point. Note that this opera- 
tion is not the spatial inversion (since all vibrations considered here form even 
representations), which are invariant under inversion, but a twofold rotation. 
Finally, for n = 5 S — |, the electronic state is orbitally nondegenerate, thus 
no JT distortion takes place. 

Table [S] collects some quantitative information about the contribution of 
each mode to the amount of JT distortion at each minimum. As expected, 
the largest distortion involves always the lowest H g mode, which is the most 
strongly coupled one (see Table HJ). Also, the D 5 distortions receive no contri- 
bution of the G g modes, which contribute to all the lower-symmetry minima 
instead. 

Table 0] contains some information about the connectivity of the minima in 
Q space. In some cases, the specification of the number of first, second, etc. 
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Table 5: The JT distortion at the minima, for each mode and value of the 
charge n and spin S. The distortions \Qi\\ are given in units of the length 
scale Xq(oj) = ^Ti/iumo) associated to each harmonic oscillator (mc is the 
mass of the C atom). The xo^a) for the G g and H g modes of Ceo are: 76.3, 
70.4, 60.3, 50.3, 46.1, 43.0; 103.7, 80.9, 62.6, 59.8, 50.1, 47.0, 43.9, 42.1 pm, 
respectively. 



neighbors of a given minimum is sufficient to clarify completely the topology 
of the minima in the 64-dimensional space. In particular, the wells of both 
the n = 2 S = and the n = 4 S = 2 surfaces are located on the six vertices of 
the five-dimensional regular simplex, the generalization of a tetrahedron, each 
minimum being equally distant from all the others: this is analogous to the 
previously determined minima of V ad for n — 1 S — | [I]. In analogy, the 
connectivity of the 10 D^d minima for n = 4 S = is the same of the one 
depicted in Fig. lb of Ref . jl] . 

For the other cases of lower symmetry, the number of neighbors of any given 
order must be complemented by some extra connectivity information. First, we 
observe that the minima for n = 2 S = 1 and for n = 3 S = | are exactly the 
same. Indeed, these two case are related by a particle-hole symmetry applied 
only to one spin flavor. For all nonequivalent cases, the complete topological 
information about the wells is contained in the connectivity matrix C(n,S), 
whose matrix elements indicate that minima i and j are C (n, 5) i -th neighbors. 
We report those matrices in the Appendix. Careful exam of C (n, S) for n = 2 

5 = 1 and n = 3 S — | shows that each of the 15 minima is linked to four 
nearest neighbor minima, which, in turn, are linked to more minima, forming 
a completely connected regular polytope. The matrices C (3, |) and C (4, 1), 
show that, for n = 3 S = | and n — 4 S — 1, the 30 minima are divided into 

6 pentagonal "clusters" of five nearest-neighboring minima. In contrast, for 
n = 5 S = |, nearest-neighbor wells come in pairs. Finally, the 30 minima for 
n = 5 S = |, show the largest connectivity, and sit at the vertices of a highly 
symmetric polytope. 

4 Vertical excitation energies 

In Table El we report the range of "vertical" excitation energies AE for all final 
spin symmetries S", in the frozen minimum configurations Q m in(^, S), for all 
values of n and S. The complete spectrum (available upon request from the 
authors) is very dense and not much informative. The listed energies give a 
quantitative prevision of the spectral range where a fast (optical) spectroscopy 
is likely to locate the intra-band HOMO excitations of the Cg^ ions. For the 
experimentally most accessible case n = 2 S = 1, here follows the complete 
list of the triplet-triplet excitation energies: 127, 149, 150, 178, 182, 218, 326, 
337, and 346 meV. 
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Table 6: The lowest and highest vertical excitation energies (in meV) calcu- 
lated assuming that the CgQ 1 " ion remains frozen in one of the adiabatic minima 
when the electronic state is excited. The first two columns fix the relevant dis- 
tortion. The third column indicates the spin S' of the excited states considered. 
The excitation energies in the last two columns are referred to the adiabatic 
energy of each specific minimum, reported in Table El 



5 Non-adiabatic corrections 



The leading quantum correction to the static JT energetics is given by the 
zero-point energy gain due to the softening of the vibrational frequency at 
the JT-distorted minima [Sj. To obtain this information, by finite differences 
we compute the Hessian matrix of the second-order derivatives of the lowest 
adiabatic potential sheet, at one of the static JT minima Q m i n 
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The vibrational frequencies u)j at the distorted point are the eigenvalues of the 
Hessian matrix. There is no square root involved in the ouj, as the coordinates 
are scaled with the harmonic length scale Xq(uj) defined in Table In the 
harmonic approximation, these "new" normal mode frequencies Uj contribute 
a zero-point energy J2j to that minimum configuration: the difference 
between this and the original zero-point energy J2iAn \Tu^iA gives the leading 
quantum correction 
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to the classical energy. We rewrite this correction as 



S) = - ^ \7~L{iAfj,}{iAn} ( n > S) — fiUJiA 
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using the invariance of the trace under change of basis. 

Table [7| displays the lowest adiabatic energies for given charge and spin 
in various approximations. The first column reports the adiabatic energy, 
-^ciass = ^ ad (Qmin), as in Table El The zero-point energy corrections in the fol- 
lowing column are comparable in magnitude to the leading adiabatic energies. 
In particular the very large values of zero-point energy gain for n = 5 S = \ is 
associated to very shallow minima, connected by low barriers. The lowest vi- 
brational frequency is as small as uj\ w 2.3 meV. In the rather close competition 
between the Coulomb physics (Hund's rules) and the JT physics (anti-Hund 
behavior) the zero-point correction is very important, and, as shown by the last 
column of Table reduces drastically the spin-gap (n = 3), or even changes 
the ground-state symmetry in favor of an intermediate in = 4) or low (n = 5) 
spin state. 
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Table 7: The total adiabatic energy -E c i ass = ^ ad (Q m in) (in meV) of the lowest 
electronic state for each n and 5*, including the e-v and e-e contributions from 
H YlD + H c -v + -f^c-c (but excluding the [Un(n — l)/2] term), for Cgg". The 
following column contains the leading non-adiabatic correction E zero , the zero- 
point energy defined in Eq. (fT4*|) . The last column reports the adiabatic energy 
-E-ciass corrected by the zero-point term E zeTO : for n > 4 it leads to a different 
ordering of the spin states. 



The zero-point correction treated here represents the g® term of a large- 
coupling expansion, where the adiabatic energy -E c iass is the leading (gf) term. 
The next corrections to be considered, of order g^ 2 , are associated to tunneling 
among minima, possibly affected by Berry phases |3J El E2] Tunneling is 
likely to be especially important between the pairs of neighboring minima of 
the n = 5 S = | adiabatic surface. Tunneling will be dealt with in future 
work. 

6 Discussion and Conclusions 

In the present calculation both e-e and e-v interactions are included for the 
HOMO shell of Ceo- E-e exchange terms are treated essentially exactly, in 
the assumptions that (i) inter-band couplings can be neglected, and only act 
as a renormalization of the Coulomb parameters and that (ii) the latter are 
independent of the charge n in the HOMO. In principle, due to both orbital 
and geometrical relaxation, the effective Coulomb interaction © will depend 



on the instantaneous charge state of the fullerene ion. However, this effect, a 
very important one in single-atom calculations, is expected to be fairly small 
in such a large molecule as C 60 . In a JT system, the coupled phonons should 
in principle be treated fully quantum mechanically, as nonadiabatic effects 
may be important. However for strong electron-phonon coupling, the leading 
terms (of order g 2 ) are obtained in the adiabatic approximation, by studying 
the minima of the lowest adiabatic potential surface V ad (Q,). Non- adiabatic 
effects are taken into account to the next order (g°) by the calculation of the 
new harmonic oscillation frequencies close to the adiabatic JT minima. These 
zero-point corrections are significantly large, and they can even reverse the 
theoretical prevision for the symmetry of the ground state of the Cg + ion for 

4 < n < 6. Tunneling matrix elements which mix different minima to suitable 
dynamical combinations restore the original icosahedral symmetry and provide 
the next-order (g~ 2 ) quantum correction to the energetics. These terms, which 
will be the subject of future work, are likely to be especially large for n = 5 

5 — |. 

The present calculation was carried out in the linear e-v approximation. As 
the coupling and thus the distortions are fairly large, quadratic and higher- 
order (in Q) couplings and vibrations anharmonicity could be important. Un- 
fortunately, no estimate for those higher-order couplings is available yet. 

The parameters used in this calculation, both for e-e and e-v interaction are 
most likely underestimated by the local density approximation used in their 
determination, as discussed in Ref . [H Ej . Consequently, both the Coulomb 
repulsion and the phonon-mediated attraction calculated within the local den- 
sity approximation are likely to need a rescaling by a similar factor of order 
two. Indeed, the balance between the two opposing interactions is delicate in 
CgJ" ions (as demonstrated by the presence of both high-spin and low-spin local 
ground states in different chemical environments [TBI EH [33 EH UH UH1 CHj)- 
In Cqq e-e interaction prevails at the adiabatic level: high-spin states are fa- 
vored, as experiments confirm for n = 2 [2.. According to our calculation, 
however, more highly charged states, close to HOMO half filling should favor 
local low-spin. 

An effective local exchange interaction favoring low spin is a crucial ingre- 
dient for superconductivity in a strongly-correlated orbitally-degenerate mate- 
rial such as a solid of doped Ceo [213 EI] • If screening and retardation effects 
could be neglected, the present single-molecule calculation suggests that super- 
conductivity should be strongly suppressed in the hole-doped solid at doping 
< n < 3, but could be recovered close to half filling 4 < n < 6. It is 
presently unclear if such a high level of hole doping is practically accessible, 
except possibly by field-induced charging [22J. 



A Appendix 



We report here the connectivity matrices whose matrix element (ij) indicates 
that minima i and j are C (n, -th neighbors, for all the nontrivial n, S cases. 
We include up to 9th neighbors, substituting those of higher order with a dash. 
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The 60 minima for n = 5 S = | are conveniently split into two sets of 30 
minima connected by the Ih operations. The minima in the "b" block are 
one by one ordinately opposite to those in the "a" block. Accordingly, the 
structure of the C-matrix is as follows: 



C 5 
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C(5,l) a C(5,|) b 
C(5,i) b C(5,l) a . 
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